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Abstract 

We investigate the quantum decoherence of frequency and polarization variables of photons via 
polarization mode dispersion in optical fibers. By observing the analogy between the propagation 
equation of the field and the Schrodinger equation, we develop a master equation under Markovian 
approximation and analytically solve for the field density matrix. We identify distinct decay be- 
haviors for the polarization and frequency variables for single-photon and two-photon states. For 
the single photon case, purity functions indicate that complete decoherence for each variable is 
possible only for infinite fiber length. For entangled two-photon states passing through separate 
fibers, entanglement associated with each variable can be completely destroyed after characteristic 
finite propagation distances. In particular, we show that frequency disentanglement is independent 
of the initial polarization status. For propagation of two photons in a common fiber, the evolution 
of a polarization singlet state is addressed. We show that while complete polarization disentangle- 
ment occurs at a finite propagation distance, frequency entanglement could survive at any finite 
distance for gaussian states. 

PACS numbers: 42.50.-p, 03.65.Yz, 03.65.Ud 
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I. INTRODUCTION 



Polarization mode dispersion (PMD) in optical fibers limits fiber performance when de- 
signing optical channels with high bit rate Physically, the origin of PMD is optical 
birefringence caused by the asymmetry of the fiber due to factors such as external me- 
chanical stress and temperature fluctuations, resulting in different group velocities for two 
orthogonal polarization modes. In addition, stochastic optical birefringence inside a single 
mode fiber leads to the random coupling of the two polarization modes, thus causing effects 
such as pulse widening , and the fluctuations of arrival times of pulses grow with the 
square root of the propagation distance [3]. Existing major applications of quantum com- 
munication, including quantum cryptography [4] and quantum teleportation , rely on 
quantum entanglement between photons as a crucial element. Therefore strategies to cope 
with decoherence have been investigated recently, for example, protection schemes based on 



decoherence-free subspace (DFS) have been proposed 



3, y, q, 



10] . Typically, a decoherence 



free two-photon state involves polarization singlet states with both photons having the same 



frequency 



llj . However, it is known that for photons with different frequencies, decoherence 



by PMD cannot be avoided due to the breaking of symmetry of the collective states 



121. 



Entangled photon pairs produced by spontaneous down- conversion can be hyp erent angled 



m 



14j, i.e. entangled in various degrees of freedom (DO 



polarization 15l. Il6j|. frequency 17|, angular momentum 



?s) of the photons. For example, 
la . Q, and energy time {20] are 



variables that can be exploited. It is thus important to address how quantitatively hyper- 
entangled photons disentangle for each DOF, by investigating the decoherence of each DOF 
separately. In this paper, we focus on polarization and frequency hyperentangled photons 
and explore how the entanglement of these two DOFs may affect each other inside optical 
fibers with stochastic PMD. 

The main purpose of this paper is to determine disentanglement length scales associated 
with frequency and polarization variables due to PMD decoherence, and to quantify the 
residual entanglement as photons propagate. To this end we will investigate three situations 
of photon propagation as depicted in Fig. 1. First we examine the propagation of single 
photon states, then hyperentangled two-photon states, in both polarization and frequency 
DOFs. In one case we consider each photon passes through separate fibers, and in another 
case both photons pass through the same fiber, all of length L. For each of these cases, we 
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22 1, which is based 



approach the problem by employing the master equation technique 
on an analogy between evolution of a state experiencing PMD along propagation direction, 
and that of spin-half particles in a stochastic magnetic field. By taking into consideration 
the frequency-dependent coupling strengths of the photon states with fiber birefringence, we 
analytically solve for the output density matrices. 

The organization of this paper is as follows. After providing the basic equations of 
our PMD model and quantum states of photons in Sec. II and III, we investigate the 
propagation of a single-photon wave packet in Sec. IV. By obtaining the single-photon's 
purity functions for each DOF, we characterize the loss of purity of each DOF of the photon 
as it propagates. The results importantly provide the characteristic decoherence lengths for 
individual photons. In particular, we obtain the pulse width at the output, which defines the 
minimum separation of well-resolved input pulses. In section V and VI, we investigate the 
dynamics of frequency and polarization disentanglement corresponding to the two-photon 
cases shown in Fig. lb and lc. Our analysis is based on Peres and Horodecki's powerful 
criterion of entanglement, known as the PPT (positive partial transposition) criterion 23j, 



24 



251 ] . and quantum entanglement is quantified by the negativity of partially transposed 



density matrix [26]. For separate fiber propagation (Sec. V), we show that frequency 
disentanglement is independent of the initial polarization status. In addition, finite length 
disentanglement is possible for both DOFs, each having distinct characteristic length scales. 
In Section VI, we address the disentanglement process in common fiber propagation, and 
particularly, we examine a polarization entangled two-photon state in the singlet form with 
Gaussian distribution of frequencies. The results provide insights about the disentanglement 
processes for states near DFS. Section VII is devoted to our conclusions. 



II. A MODEL FOR STOCHASTIC POLARIZATION MODE DISPERSION 

In this section, we present basic equations for an optical pulse in optical fibers subjected 
to stochastic PMD [27j . In a single mode fiber, only two modes with orthogonal polarizations 
are supported. For an optical pulse propagating along z in a linear and birefringent medium, 
fiber birefringence b(a>, z) = 6161 + 6262 + 6363, where ej are unit vectors in the Stoke's space, 
alters the polarization the field depending on its frequencies. The birefringence effect on the 
field modes, each of frequency u, can be derived directly from the wave equation, which 
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a) 



b) 



c) 




|0> |0> 



FIG. 1: Single photon pulses passing through fibers of length L along z direction: a) A single- 
photon input pulse; b) two-photon entangled input pulses each passing through separate fibers; 
and c) two-photon entangled input pulses passing through a common fiber. 

gives 

-^E(cu,z) + k 2 (u J ,z)B(u J ,z) = 0. (1) 

Here k 2 (u, z) = [3 2 (uj, z) + /?o(^, z)b(u, z) ■ er is the propagation tensor, with (3q{u, z) as the 
common propagation constant, and the vector a is formed by the Pauli matrices given by 
a = <xiei + <7 2 e 2 + (T-ie 3 , with 




Extracting the fast oscillating common phase from the field, we define E(c<j, z) = 
E(u)exp [— i J Q Z Po(uj, z')dz'\ A(uj,z) with E(u) as the amplitude of the mode and A(u,z) 
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as its frequency spectrum of the 2-dimensional Jones vector A(t, z) 



A(o>, z) = I dtA(t, z) exp(iu>t). (3) 



'2tt 

We adopt the adiabatic approximation, assuming that the birefringence vector b(u;, z) and 
thus the polarization of the pulse varies slowly along z PJ. As a consequence the term 
Jp-A(u;, z) becomes negligible, and the effects of birefringence on the spectrum A(u, z), as 
the pulse passes through the fiber, can be described by the following first order differential 
equation [l| 

d 1 

i—A(u,z) = -[b(co,z)-(T}A(cu,z). (4) 

Note that for non-dispersive channels, b(uj,z) ■ a gives eigenvalues ±uj[riF(z) — ns(z)]/c 
determined by the refractive indices associated with the fast mode np and slow mode rig- 

It is important to note that Eq. (j4]) is a kind of Schrodinger equation with the 'Hamilto- 
nian' h — | [b(uj, z) ■ <r], in the form similar to that of a spin-half system interacting with a 
magnetic field. For deterministic evolution, we can express the output pulse with a unitary 
transformation, 

A(u,L) = e- if ti hdz A(uj^) (5) 

with T referring to a position ordered integration analogous to the time ordered integration 
in quantum theory. In terms of column vectors in Jones space, we have 

A(w,0)=| I and A(u,L) = I I, (6) 




where C™ and C° ut (J — 1, 0) are polarization amplitudes for input and output fields obeying 
the normalization condition: \C™\ 2 + |C} n | 2 = | Cq ut 1 2 + I C° ut | 2 = 1. For conceptual clarity, 
we will assume that the regions z < and z > L are free of birefringence. In this way, 
A(lu, z < 0) = A(co>, 0) and A (a;, z > L) = A(u, L). 

The stochastic nature of our PMD model originates from the randomness of h(u, z). In 
this paper we adopt the assumption from [2i], that the randomness is due to the fluctuations 
of optical axis along the fiber. In addition, the frequency dependence of b(ui, z) is determined 
by fiber material properties only and therefore should not change over different positions. 
In this way it is plausible to assume ^]: 

b(u,z) = f(u)b(z) (7) 
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where h(z) is stochastic, and the f(u) is a deterministic function defined by the material 
and it can be expressed in terms of a Taylor series about the peak frequency ujq of input 
pulses: 

f{uj) = ju + suj(uj — Wq) H • (8) 

Note that in non-dispersive media, only the first term remains and most of our discussions 
below will be based on such an approximation. To specify the statistics of h(z), we assume 
;hat h(z) is a random process with zero mean, and it has the two-point correlation function 

b(zi)b(z 2 ) = r ] 2 Yd(z 2 - z x ). (9) 

Here the bar refers to the ensemble average, and rj characterizes the average strength of 
birefringence. We remark that the assumption of delta correlation function ([91) is not a 
strict requirement. As long as the correlation length of h(z) is sufficiently short, in the sense 
that A(u>, z) does not change significantly within the correlation length, then it is justified 
to employ the Markovian approximation in obtaining the master equation in later sections. 



III. GENERAL DESCRIPTION OF SINGLE AND TWO-PHOTON STATES 

Let us first examine a deterministic situation corresponding to a given realization of 
b(u,z). In the previous section, we have seen that if an incoming wave of frequency to 
incident from the left with the polarization state (Cf ,Co n ), then according to Eq. (jHJ) and 
Eq. (El), there is an outgoing wave propagating to the right with the polarization state 
(C° ut ,Co ut ). Since the system (or the wave equation) is linear, an incoming single-photon 
follows the same transformation rule to become an outgoing photon. Let 4>(u) be the fre- 
quency envelope of the input single-photon wave packet, then the input and output state 
vectors, denoted by \^^) and l^out)* take ^ ne f° rm: 

I**?) = J du<f>(u) |w> ® [Cj» |1> + C » |0>] (10) 

|*St> = / du<t>(u) \u>) ® [C^iu) |1) + |0>] (11) 

where \u>) is the frequency basis vector defined in the birefringence free [h(u, z) = 0] system, 
and |1) and |0) respectively correspond to horizontal and vertical polarization basis vectors. 
We point out that in writing Eq. (10) and (11), we have employed a rotating frame such 
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that the phase factor e~ tuJt due to the free field evolution of \u) has been removed. This 
is equivalent to the representation in interaction picture. If Schrodinger picture is needed, 
we just need to replace \u) by e~ luJtin \uj) in Eq. (10), and \u) by e~ luJtout \uj) in Eq. (11), 
with t- in and £ out are instant of times defining the input and output states. Both £ in and £ out 
should be chosen in such a way that the input and output wave packets are far away from 
the birefringence interaction region. 

It is important to note that if we treat the interaction length L = z (Fig. 1) as a 
parameter, and let l^outO 2 )) be the output state corresponding to a birefringence fiber of 
length z, then l^outO 2 )) * s governed by the Schrodinger-like equation according to Eq. (TJ]), 
i.e., 

^\^)) = n^(z)\^ t (z)) (12) 

where z plays the role of time, and 

H {1 \z)= j ' du\u) (u\®^[b{u,z) ■ a] (13) 

plays the role of Hamiltonian. 

In the case of two-photon states we will restrict our discussion to systems involving 
two distinct single-photon pulses, A and B, such that each pulse contains a single photon 
(Fig. 1). In other words, we can label the photons as two subsystems A and B. The 
distinguishability of the two photons can be achieved in two physical situations of interest 
here. The first situation is illustrated in Fig. lb in which the two single-photon pulses 
individually propagate in two different optical fibers, and the second situation is when two 
spatially (or temporally) separated photons propagate in the same fiber (Fig. lc). In both 
cases, we have the input-output state vectors: 



sa,sb=0,1 

with <J)(ua, Ub) being the normalized frequency envelope of the input 2-photon wave packet, 
i.e., J J dujAduJB^i^A-, ^s)| 2 — 1; an d the (u) and C°"* s describe the joint polarization 
amplitudes for input and output states at the corresponding frequencies. We remark that 
input states with non-factorizable 4>{uj Ai ujb) correspond to frequency entangled states, and 
similarly, non-separable C™ ASB means polarization entanglement. 



= / du A du B <j){uJ A ,ujB)\uJ A ,ujB) ® ^2 C T a s b ( uj a,^b)\sa,Sb) (14) 

"* s A ,s B =0,l 

l*£t> = / du A du B (f>(uJA,UB)\uA,UB) ® C sTsb( uj a,u b )\sa,sb) (15) 



The fact that the two distinguishable single-photons do not interact allows us to treat their 
evolution by the transformation rule as in the case of single photon. Similar to Eq. ffl2|) . we 
can treat the interaction length L = z as a parameter and obtain the Schrodinger equation, 

i§-yi(z))=rt 2 \z)\^i(z)} (16) 

with 

H^(z) = \J du [bi(w, z) ■ trW ® \u) A (oj\ + h 2 (uj, z) ■ tr™ ® \u) B (u\] (17) 

for the case of separate fiber (Fig. lb), and hi(cu,z) and b2(oo,z) are birefringence vectors 
of the two fibers. The cr^ and cr^ are Pauli vectors for photons A and B. 
In the case of common fiber (Fig. lc), we have 

rt 2 \z) = X - J dub(u;, z) ■ [a™ <g> \u) A (u\ + * {B) ® \u) B . (18) 

which indicates that both photons experience the same birefringence interaction. Note 
that the pure state vectors discussed above can be considered as quantum trajectories, 
corresponding to a single realization of a birefringence vector. To address the stochastic 
problem, we will need to to perform averaging via the density matrices. In later sections, 
theoretical analysis of the density matrices will be carried out according to the master 
equations for each of the cases in Fig. 1. 



IV. DECOHERENCE OF A SINGLE-PHOTON STATE 



Let us first discuss the master equation describing a single-photon state passing through 
the fiber. From Eq. (fl2l) . the density operator ^(z) of the state obeys the equation 



dz 



(19) 



where the subscript s stands for a given realization of birefringence defining 7i,^(z), i.e., 
before ta king the ensemble average. We follow the standard strategy to obtain the master 
equation [221 ]. 

.2 3 



d_ 

dz 



(i) 
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i=l 



>ri - f,f ip {1) - £ (1) r, : r 



p 



(20) 
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where we have used = pi in order to simplify the notation, and f j = J du>f(u) \u>) (u\® 
&i is defined. For later purposes the single-photon density matrix can be expressed explicitly 
by, 




dujdui' p 8a i (u,u';z)\w,s) {w',s'\ (21) 



■;S' = 0,1 



where p ss i (u, u/; L) are matrix elements. 

We point out that the derivation of the master equation Eq. (j20p is based on the Bloch- 
Redfield-Wangsness approach known in nuclear magnetic resonance literature [221 ] . Alter- 
natively, the same master equation can be derived by treating the fiber medium as a bath 



with many degrees of freedom 12j, |2lj. The former approach, which we adopt here, can be 
understood more transparently by considering the fiber as composed of concatenating un- 
correlated short sections of length Az, and Az is set to be long compared with the coherence 
length of h(z), but is small so that the change of state can be approximated by keeping the 
Dyson series up to the second order in Ti,^ . Then by Markovian approximation Eq. 
the master equation ( 1201) in fact corresponds to the 'coarse rate of variation' Apy/Az upon 
ensemble average. 

Now, we consider a general single-photon pulse which is initially polarized along |1), 

|* in ) = |1) <g> J duo(j){uj) \u) . (22) 

We remark that the result is the same for arbitrary polarization direction due to the sym- 
metry caused by the randomizing effect of the birefringence fluctuations, and hence we set 
it to |1) for convenience. The input pulse envelope is set as a Gaussian wave packet 

1/4 M^o) 2 " 



M = ( ^ ) ex P 



K 2 



(23) 



where k indicates the width of the Gaussian envelope and the peak frequency ujq is in the 
optical range. With the input condition p in = |1) (1|<8>/ / dujduj'(j){uj)(j)*{uj') \uj) (u'\, we solve 
the output state governed by the master equation Eq. (TSUI) , and the details are presented in 
Appendix (jA]). The solution for the density matrix is given by, 

Pll (u,u';L) = V) K AlL + e~ A2i ) 

p 00 (^, cu'; L) = ^(u;)0V) (e~ AlL - e~ x * L ) , (24) 



where the values of Aj > at given to, uo 1 are 

= ^[/H-/(^)] 2 

A 2 (a,,u/) = ^[3/H 2 + 3/(^) 2 + 2/M/(^)], (25) 

It is interesting to note that Ai is a difference of the frequency profiles and A2 is a sum, giving 
Ai <C A 2 in the optical region. The off diagonal elements are pio(ui, u?\ L) = Pq\{u, uj'\ L) = 0. 
Having solved the matrix elements, the output density matrix p^'(L) can therefore be 
obtained according to Eq. (|21[) . It can be observed that in the long length limit L — > 00, a 
complete depolarization occurs, since only diagonal elements pn(u,u;L) = p 00 (tj,u>; L) = 
\ I 4>{uj) 1 2 remain. 



A. Pulse spreading 

Let us introduce the quantized field operator, 

£{z,t) = J duo e^ (z - ct)/c a w + h.c. (26) 

where is the corresponding annihilation operator. To visualize the change of output pulse 
shape due to PMD, we examine the intensity operator I = S^S by evaluating its expectation 
value with respect to the output pulse. The calculation is quite tedious, and we present the 
results only. For linear non-dispersive fibers f(oo) = ju, we find that the spatial dependence 
of the averaged output intensity associated with the two polarization states are given by 

8L/L C 



<jl I g 2[\+6L/(L c J I y\ g 2[\+2L/(L c v 2 )} l+iL/(L c v z ) 



2 I v/l + 6L/(L c i/2) v /[l + 2L/(L cI /2)][l + 4L/(L c z/2)] 

/{ _ K 2 T 2 _ K 2 T 2 _ SL/Lc 

g 2[l+6L/(L c i/ 2 )] g 2{l+2L/(L c v 2 )] l+4L/(L c v 2 ) 

y/l + 6L/(L e v*) ~ y/[l + 2L/(L e v*)][l + 4L/(L e u*)] ? ^ 

where r = z/c — t, with t being a sufficiently long time so that the entire pulse has exited 
from the birefringence fiber of length L. In addition we have defined the characteristic 
decoherence length 

w^y 1 , (28) 



and the dimensionless ratio v = Uq/k. 
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In the narrow bandwidth case where v 3> 1, the second term in Eq. (1271) decay approx- 
imately exponentially with the decay length L c /8. Such a decaying length scale is much 
shorter than that in the first term. Therefore as L ^> L c increases, only the first term 
remains, equalizing both elements and hence showing depolarization. It can be seen that 
the width of the output pulse is approximately Ca/2[1 + QL/ (L c is 2 )}/k, which has a VZ de- 
pendence when L ^> L c v 2 is sufficiently large. Such a vT dependence were also reported 
in general classical consideration j^]. The pulse width at the output is greater than that of 
input, thus sets a lower bound to the distance between input pulses in order to allow the 
output pulses to be non-overlapping, or distinguishable. 



B. Purity functions 

Next, we investigate the decoherence of the single-photon pulse through its purity defined 
as fi = Tr(p 2 ) < 1. Purity is closely related to the reciprocal number of effective modes 
required to contain the whole state, hence the equality sign holds only when the state is 
pure, i.e. can be completely represented by one mode of the basis. We calculate the frequency 
purity, the polarization purity and the overall purity. 

To obtain the frequency purity, we trace the polarization freedom of (L) and consider 
only the frequency DOF. From Eq. (j2li) and Eq. ( |24|) . noting that 

Tr s (p (1) (L)) = J J dujduj'(t){u)(t)*{uj')e- XlL \u) (u'\ , (29) 

the frequency purity is given by, 

^(L) = Tr{[Tr s (p«(L))] 2 } = =L= (30) 

a/1 + QL/(L c is 2 ) 

which decays algebraically with L. Note that fM u (L) depends on a decay length scale z/ 2 L c , 
meaning that ^(L) has a slower decay rate for a greater u, i.e., narrower spectrum. A 
remark is that the effective number of frequency modes (as measured by p" 1 ) required to 
represent the state increases from 1 to infinity as L increases. In fact, /i" 1 shares a similar 
functional form with the output pulse width described in the previous section. This suggests 
that the pulse widening can be a measure of PMD decoherence of frequency variables for 
Gaussian initial states. 
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Next we find the polarization purity of p^(L) after tracing the frequency freedom from 



p {1) (L). Noting that 

Tr u (pF>{L)) = J du[ Pll {u,u;L)\l)(l\+p O0 {u,u;L)\0} 
the polarization purity p s {L) is therefore 

p s (L) = Tr{[Tr w (pW(L))] 2 } 



(31) 



dujpw (cu, u; L) 
1 



+ 



dujpoo (uj,uj; L) 

l(iL/Lc 



1 + 



_g 1+4L/(L C ,^) 



i + AL/{L^y j' (32) 

decreasing from 1 to the minimum value 1/2, meaning that the state is more spread out in 
the 2-dimensional Jones space to the fully mixed situation as fiber length increases. Note 
that the decay length scale for p s {L) is L C /1Q in the narrow bandwidth case with v ^> 1, 
which is shorter than that for p^iV). In addition, appearance of the exponential factor in 
Eq. (13"2"|) indicates a faster decay rate than that of p^L). 
The total purity is also found by 

pitotai(L) = Tr{[p«(L)] 2 } 

1 



16L/L C 



g l+iL/iLcv 1 ) 



2 I ^l + 6L/(L c v*) v /[l + 2L/(L c z/ 2 )][l+4L/(L c z/ 2 )] 



(33) 



again having a slower rate of decay for a greater z/, and obeying p^L) ■ p s (L) < p tota i(L). 
We conclude that decoherence for single-photon state in Eq. (|22|) is complete only upon 
L — > oo. 



V. DISENTANGLEMENT OF TWO-PHOTON STATES IN SEPARATE FIBERS 

In this section we study the disentanglement of an entangled two-photon state propagat- 
ing along separate fibers, each with random birefringence bi(z) for (i = 1,2). We assume 
that they are made from the same material, hence obeying (bj(zi)bj(z2)) = r] 2 I 5ij5(z2—zi). 
The 'Hamiltonian' in this case is given by Eq. (fT7|) . With similar assumptions as in the 
single-photon situation, we obtain the master equation governing the ensemble averaged 
density matrix p^ = pf^ 

^ (2) = t £ E [af ^f?> - f ? } f ! j) p (2) - P (2) f ! j) ff 1 , (34) 



dz r 4 

j=A,B i=l 
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where = J duf(u) \uo)- (uj\ ® o" f for the two photons j = A, B. It can be noted that 
the master equation consists of two decoupled parts, each for an individual fiber. As in the 
previous section, it will be convenient to write the two-photon density matrix explicitly 



pP\z)= / / duj A duj B duj' A duj' B ^ Ps A s B s' A s' B (wa, ub, u' A ,d B \ z) 

x \ua, uj b ) (u' a , u' B \ <S> \sa, s B ) (s' a , s' b \ . (35) 

with the matrix elements p SA s B s' A s' B (&A, wg, u' A , uj' b \ z). 

In this section we investigate an initially hyperentangled state a polarization Bell 
state with frequency entanglement: 

\Mn) = J J <fi(uj A ,ujB)duj A dujB \ua,Ub) ® |V>Bell) , (36) 

where the frequency envelope is assumed to be a double Gaussian, with a peak frequency 
Uq, as follows: 



<fi(uj A ,uj B ) = W^^exp[-a 2 (^A - wb) 2 - (3 2 (uj a + u B - 2o> ) 2 ], (37) 
I if 

where the width a > (3 corresponds to a more frequency correlated state and (3 > a indicates 
a more frequency anticorrelated state. We consider an input singlet pulse I^Beii) — (| 10) — 
1 01)). The evolution of the four Bell states, including the singlet state and the triplet states, 
follow same calculation steps and thus only the singlet state evolution is discussed in detail, 
presented in Appendix (jB|) . 

Following the master equation Eq. fl34|) . the evolution of the only non-zero density matrix 
elements at u A , ljb, uj' a , oj' b can be found as 

Pmi = Poooo = ^t4>(uja,ub)4>*(ua,uj b ) (e' ClL - e~ UL ) 

Pioio = Poioi = j<P{^a, u B )4>*{u' A ,u' B ) (e _ClL + e" C4i ) 

Piooi = Pono = -^4>{^A,^B)<P*{uj' A ,uj' B )e~' :4L , (38) 

where the values of Q at given u A , ub, ui' a and u' B are shown in Eq. (1B4[) . We can see 
that (i <C ^4 i n the optical region. We remark that the Bell states have zero projection 
to the spaces characterized by £2 and (3 and therefore (2 and £3 do not contribute to the 
state evolution, which is shown in Appendix ([B]). In particular, d = only for the diagonal 
elements, and thus the steady state is a completely depolarized one, with p m i = pioio = 
Pom = Poooo = i \4>{ua,u b ) I 2 - 
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A. Characterization of entanglement in terms of negativity 



24 



According to Peres and Horodecki's PPT (positive partial transposition) criterion 23, 



25], if the partial transposition of a bipartite density matrix (denoted by p Tjv ) has one 



or more negative eigenvalues, then the state is an entangled state. The negativeness of p TA 
turns out to be a necessary and sufficient condition of mixed state entanglement for two- 
qubit states and bipartite Gaussian states. To quantify how much entanglement survives 
by PMD decoherence in polarization and frequency variables, we calculate the negativity of 
the corresponding DOFs. 

The negativity Af of the state p is defined by [26 ] 

1 1^1 1 - 1 

Af = ^—^ — - (39) 



which is an entanglement monotone under local operation and classical communication 26]. 
Specifically, the polarization negativity is found by first tracing the frequency variables and 
taking the trace norm of the partial transposition (sa *-> s' A ), i.e. 

u, = ( 40) 

Similarly, the frequency negativity is obtained by finding the trace norm of the partially 
transposed density matrix with the polarization variables traced, i.e., 

K = \\{^,,Ap<W*\\-\ (41) 

Equivalently, negativities can be obtained by summing the absolute value of negative eigen- 
values of the partially transposed matrices 28] , which we adopt in the following sections. 
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B. Polarization disentanglement 



Now we discuss polarization disentanglement of the two-photon output state by finding 
the negativity M a of the state. We first obtain 



{^ a ,.Ap {2 \l)]} Ta 




duJAdtxlB 



( Pun \ 

Pioio Piooi 

Pono Poioi 

V poooo/ 

(kiX-x) o o 

o \{i + x) o 

i(l + x) o 

V -\x o o \{i-x)J 



(42) 



where by assuming f(uS) = 70;, 



X = exp 



-16L/L C 



l + {2L/(3 2 ulL c )_ 
The polarization negativity M s can thus be found as 

3 




(43) 



4 



X 



1 

4' 



(44) 



Note that Eq. (1441 is applicable to all four Bell states since we have the freedom to redefine 
the polarization bases and phases in the second fiber due to the symmetry caused by the 
stochastic birefringence fluctuations. We also remark that M s has a decay length scale 
L c /16, and we note that the same decay length scale exists for polarization purity in the 
single photon case. In addition, we observe that finite length disentanglement is possible 
when ^ < 1/3, which can be solved numerically. As the exponential factor of x is depending 
on the value of (3, we fix a and plot the trend of polarization negativity with varying (3uq 
in Fig. [21 We see that in general the critical disentanglement length increases and peaks at 
some value of /3ujo, then level off to some finite value as j3uo tends to infinity. 
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FIG. 2: The polarization negativity for the two-photon separate fiber case for aujQ = 1000 with 
varying /3 and L. The black solid curve indicates the critical disentanglement length. 

C. Frequency disentanglement 

To calculate the frequency negativity, we trace the polarization variables as follows: 



It is insightful to note that this equation holds not only for Bell states, but for any general 
initial polarization states having the form of Eq. (1361) {^J. An important consequence is 
that frequency entanglement in separate fibers is independent to initial polarization status. 
This result is in contrast to the polarization negativity we presented previously, which is 
dependent on the initial frequency envelope shape. 

The partial transposition of the frequency density matrix can be found by 
P Ta (uja, ujb, uj' a , uj' b ) = (u' a ,ub\^s a ,s b \p^ 2 \L)]\ua,u' b ). For linear non- dispersive medium 
f{uj) = ju, the frequency density matrix remains Gaussian. This allows us to determine its 
negativity can be analytically from the general formula for Gaussian states in 30] , 




(45) 




for a 2 ulL c > (3 2 uj%L c + 3L 
for !3 2 ulL c > a 2 ulL c + 3L 



(46) 



It is interesting to see that finite length frequency disentanglement occurs when 



, , 2 

L > —\a - (3 \L, 



(47) 
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which is universal to any initial polarization state. The expression signifies that the system 
is 'less robust', i.e. having the critical length of disentanglement tends to zero, if the initial 
frequency envelope has a ~ (5. Note that the polarization entanglement mainly depends on 
the value of /?, while the frequency entanglement increases with the difference between a 2 
and (3 2 . 



VI. DISENTANGLEMENT OF TWO-PHOTON STATES IN A COMMON FIBER 



Next, we examine the decoherence problem for two photons propagating in the same fiber 
of length L. Physically, we can distinguish the two photons by spatially separate them by a 
small distance, so small that the two photons still experience the same stochastic interactions 
with the fiber birefringence, yet far apart enough to be distinguishable, depending on the 
width of the output pulse which is discussed quantitatively in the single-photon section. The 
'Hamiltonian' in this case is given by Eq. ffl8l) . The master equation is as follows: 



1 = t ' 

dz P 4 

i=i 



6 

[2A,P (2) A, - kkf> {2) ~ P (2) kk] , (48) 



where 




A; = / / du A du B [f{u A )a\ A + f{u B )a\ B) ] ® \u A ,u B ) (u A ,u B \ . (49) 



We remark that the form of this equation is similar to the single photon master equation 
Eq. ( 1201) . but with the collective operator Aj. 

It is known that two photons of same frequencies in singlet polarization states lie in 



DFS, and their resistance to decoherence has been verified experimentally llj. The reason 



behind the decoherence free effect is that they are eigenstates of iiS 2 ^ with an zero eigenvalue, 
meaning that these state, do not evolve as they propagate. The hey is the identical conphng 
with the environment for both of the photons [12j. In practice, however, photons generally 
have fluctuations in frequencies variables. In the case of down conversion systems, the 
two photons are in fact anti-correlated due to energy conservation. It is thus of interest 
to examine the robustness of entanglement for singlet polarization states with a general 
Gaussian frequency envelope (f>(ui A ,uj B ), 



= / / <\>{yj A ^ B )e- WAZ « tc duj A duj B \uj A ^ B ) 



^(|10>-|01» 



(50) 
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where 4>{ijJa, &b) is a double Gaussian given in Eq. (j57|) . and the phase factor e~ tu)AZ °/ c is 
added in order to displace the peak position of the photon wave packet A by a distance 
Zq relative to B. The choice of separation zq should be large compared with the width of 
the wave packets so that the two photons can be treated as distinct subsystems, but small 
enough so that both photons would experience the same stochastic birefringence. Note that 
the displacement of photon A is simply achieved by a local unitary transformation operator: 
exp(— iz J -\oj) a(uj\cIuj) , and it commutes with in Eq. (TT81) . Therefore the phase factor 
do not affect the entanglement. For convenience, we will simply absorb the phase factor into 
the definition of \co)a in the later calculations. 

Following the master equation in Eq. (jUJ), the singlet state evolution at to a, ojb-, w' a , w' b 
are found by the steps in Appendix (jC|) . and are presented as follows. 



where 



/ 



Pxxxi Pino P1101 Pnoo 

PlOll P1010 PlOOl PlOOO 

Pom Pono Poioi Poioo 

V Poon Pooio Poooi Poooo / 



/I 



v/3 



V2 











5^1 + 2^2 -5^1 + 273^ 



-^1 + ^2 



2 

l^+2V3 





V-2 



\ 






(u'a,u' b ) 



-V w + V 



22 



+ 1 



(1) 



11 



V, 



(1) 



22 



4K 



12 



(v {1) - V 



(i) 

22 



+ 4K 



12 



V2 



V£U{UA,UB)<f>*(UA>U > B) 



it 



(51) 



(52) 



having £j > denoted in Eq. ( 1C8I) . with £i ^> £ 2 in the optical region. 

An important situation to notice is when £ 2 = 0, where no decay occurs and the matrix 
elements can survive to infinite fiber length. It can be noted that £2 = only for three cases, 
when oja = oj' a and u B = u' B , uja = ojb and oj' a = u' B , or uja = oj' b and u' A = ojb (Appendix 
C). Thus we obtain the output state in the L — > 00 limit, with the only nonzero terms as 
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follows: 



[ Pun Pioio Poioi Poooo Piooi Pono ) T 

= i0(cJ A ,CU B )0*(^,^ B ) ( 1 1 1 1 0) 

( Pirn Pioio Poioi Poooo Piooi Pono ) T 
= ]^)(u)a,ua)4>*{u a ,Ua) ( 110-1-1 



L^oo 



( Pllll Pioio Poioi Poooo Piooi Pono ) 
= i0(w A ,cu B )0*(cj B ,cc; A ) (-1 o -1 -1 -l) (53) 

We remark that the first case refers to a completely depolarized situation in which the 
density matrix contains only the diagonal elements. Furthermore, the second case refers to 
the decoherence free situation with both photons having the same frequencies, and therefore 
experiencing collective decoherence. In this case it does not decay and remains as a singlet 
state for any fiber length L. 



A. Polarization disentanglement of the singlet state 

Now we discuss polarization disentanglement of the two-photon output state by finding 
the corresponding negativity M of the state. We first obtain, 



{Tr^ B [P (2 \L)]} TA = 



where by assuming f{uj) = jiv, 



( \{l-v) -\v \ 

\H + v) 

\{} + v) 

V -¥ \{l-v)) 



v = j J du A du B \(p(uJA,uJB)\ 2 exp {-2if[f(uj A ) - f(u B )] 2 L} 



= 1 / Wl 



AL 



a 2 u$L c 



The polarization negativity is therefore 



Kf 3 1 



(54) 



(55) 



(56) 
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Finite length disentanglement occurs when 

L > 2a 2 ujlL c . (57) 

We remark that polarization disentanglement is 'more robust' if the two photons have more 
correlated frequencies, i.e., a larger value of a. In the limit a — > oo, negativity never decay, 
which is the original DFS case where both photons experience collective decoherence. In 
addition, comparing with the polarization negativity of common fiber with that of separate 
fibers as in Eq. (|44p . we see that the latter has in general a greater rate of decay due to the 
exponential decay factor in Eq. ()43l) . 



B. Frequency entanglement of the singlet state 

To investigate frequency entanglement, we first trace the polarization variables of the 
output density in Eq. (|5T|) . and take the partial transposition u A <-> uj' a , obtaining the 
density matrix elements as follows, 

pI a (u a ,uj b ,uj' a ,uj' b ) = {u' a ,ujb\^s a ,s b [p {2 \^)}\^a^' b ) =Vi + V3v2. (58) 

With a non- Gaussian form Eq. (|58|) . there is not a generally agreed analytical measure for 
mixed-state entanglement. We thus attempt to detect the presence of entanglement by 
evaluating the positivity of the partially transposed density matrix in Eq. ( J58l) . We test 
its positivity by noting that for the whole matrix to be positive, any 2x2 submatrices 
composed by extracting 4 points from two of the rows and columns of the density matrices 
should be positive. Here we will extract the points where £2 = so that they do not decay 
and are significant upon the long length limit. For the frequency correlated case we select 
uj a = u B and u' A = u' B , the second case of the long length limit in Eq. (I53I) . giving the partial 
transposed elements p^ A (u)B, uj a , u a , ujb), we consider the 2x2 submatrices for u A < ujb'- 

pl A (u a , uj b , uj a , u b ) pl A (uja, ujb, uj b , uj a ) 
pL a (^b, uj a , uj a , lo b ) pI a (uj b , uj a , lo b , uj a ) 
Checking the positivity of the matrix, we employ the fact that the submatrix is negative 
if and only if the product of off-diagonal terms is greater than the product of the diagonal 
terms. Noting that 

Pw A ( U A, uJb,uj A} (j b )Pu A ( uj b, uj a ,uj b ,uj a ) 



(59) 



pl A (uJ A , uj b , uj b , uj A )pZ A (uj b , uja, uja, ujb) 
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exp[-4(a 2 -(3 2 )(uj a -lobY}, (60) 



we see that all these submatrices are negative when a > (3, independent of the birefringence 
interaction length L. This indicates that for singlet polarization states, frequency entan- 
glement exists at any finite distance L if the frequencies of the two photons are initially 
correlated in the form of a double Gaussian. 

On the other hand, if the two photons's frequencies are initially anticorrelated, i.e., f3 > a, 
we select uja = u'b an d 0Jb = w' A for the sub matrix positivity test. This is guided by third 
case in Eq. fl53|) . We therefore consider the 2x2 submatrices 



which are negative if 



pI a (uja, ua, ua, ua) pI a (^a, ua, ub, u b ) 
Pu A {ub, Ub, u a , va) pI a (ub, Ub, u b , u b ) 



{3>^a* + g(L), 



(61) 



(62) 



where 



1 



In 



1 _|_ e 8(i^A-^s) 2 i/^o L < 



is defined. Particularly in the long length limit, Eq. 



> 0, 



(63) 



621) reduces to (3 > 



a 2 + 



4(uA-Wfl) 2 



In 4, meaning that entanglement persists at long distance if j3 is sufficiently 



large. 



VII. CONCLUSION 



To summarize, we discuss quantum disentanglement of frequency and polarization vari- 
ables, for photons propagating through fibers with stochastic PMD. Observing the analogy 
between the wave propagation inside the fiber and the Schrodinger equation in quantum 
theory, master equation method is adopted to analytically solve for the field density ma- 
trix. In this paper we investigate the single-photon and two-photon cases. For the single 
photon case, purity function for each of the DOFs is analytically calculated, quantitatively 
determining the degree of mixing, which reveals that complete decoherence is possible only 
for infinite fiber length. Pulse width of the output pulse is also evaluated, determining the 
minimum separation of pulses for them to be distinguishable at the output. Next, for en- 
tangled two-photon states with each photon propagating through a separate fiber, we show 
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that entanglement associated with frequency and polarization variables can be completely 
destroyed after distinct finite propagation length scales. Specifically, for the hyperentangled 
state Eq. (1361) . condition of polarization disentanglement is found as x < 1/3 where \ is 
defined in Eq. ( I43p . and the condition of frequency disentanglement is given by Eq. (JHJ). An 
interesting fact that frequency disentanglement does not depend on the initial polarization 
status is also revealed. For a singlet polarization state propagating through a common fiber, 
we show that polarization disentanglement in finite length is possible, though having a much 
longer critical length of disentanglement than the separate fiber case. We also consider the 
frequency entanglement in common fiber, observing its dependence on the initial frequency 
envelope. On one hand, for the frequency correlated parts of the density matrix, entangle- 
ment persists if it already exists at the input, explainable by the DFS. On the other hand, 
entanglement can manifest in the anti-correlated parts of the density matrix if initially the 
frequency envelope is sufficiently anticorrelated. 

Acknowledgments 



This work is supported by the Research Grants Council of the Hong Kong SAR, China 
(Project No. 401307). 

APPENDIX A: SOLUTION OF MASTER EQUATION FOR SINGLE-PHOTON 
STATE 

We outline the solution of the master equation Eq. ( 120]) in the following. From Eq. ( l20l . 
the matrix elements p u (uj,u'; z) and p 00 (u,uj'; z) are coupled by 




(Al) 



where 




3[/H 2 + /(^) 2 ] + 2/M/K) 

4/(w)/(w0 



4/(w)/(W) 
3[/H 2 + f{u'f\ + 2/(w)/(*/) 



) 



• (A2) 
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Diagonalizing Mi, its eigenvalues are — Aj < as defined in Eq. fl25l) . and hence we can find 
the solutions by evaluating 



Pn 
Poo 



exp(M^) 



Pn 
Poo 



(A3) 



2 = 



On the other hand, the off diagonal elements p w (u,u'; z) = Pqi(uj, u/; z) decay individually 
with the same form 



9 

—p 10 (oj,uj; z) = -X 2 p w (uj,uj'; z) 



(A4) 



Putting in the input conditions pu(uj, u/; 0) = (f)(u)(j)* (u') , poo(^,^';0) = Pio(lo, u/; 0) = 
Poi(^,^';0) = 0, we get the solutions in Eq. ffMty . 

APPENDIX B: SOLUTION OF MASTER EQUATION FOR TWO-PHOTON 
STATE IN SEPARATE FIBERS 



Following the master equation Eq. (1341) . we found that the matrix ele- 
ments p ini (u A ,u B ,u' A ,u' B ;z), p ww (uj A ,uj B ,uj' A ,uj' B ;z), p 0W i(u) A , u B , u' A ,u B ; z) and 
PooooO^a, wg, u' A , uj' b ; z) are coupled as follows: 

/ Pun \ ( Pnn \ 



d_ 

dz 



Pioio 
Poioi 
\ Poooo / 



M, 



Pioio 
Poioi 
V Poooo / 



where 



/ 



M 2 



rrii rri2 r«3 



\ 



with the elements 



777,2 mi 777,3 
77i3 mi ™2 
\ 7773 7772 TTlx J 



(Bl) 



(B2) 



mi = -3[/(^) 2 + /(^) 2 + /KJ 2 + ^^ 
777 2 = Af{uj B )f{uj' B ) 

m 3 = 4f(u A )f(u/ A )' (B3) 
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Diagonalizing M2, its eigenvalues are — Q < where 



Ci = \{mM - /K)] 2 + [3(/( Wfl ) - /K)] 2 } 

C 2 = ^{3[(/(o; A )-/(^)] 2 + [3/(a; B ) 2 + 2/(a; B )/(^) + 3/(^) 2 ]} 

Cs = ^{[3/K) 2 + 2/(^)/K) + 3/K) 2 ] + 3[(/K)-/(^)] 2 } 

C4 = ^{[3/K) 2 + 2/(^)/K) + 3/K) 2 ] 

+ [3/(^) 2 + 2f(u B )f(u' B ) + 3/K) 2 ]}, 



(B4) 



with d having the smallest value and £4 the largest. The corresponding eigenvectors are 



1 
1 

V 1 / 



^2=2 



1 

-1 

V 1 / 



^3=2 



1 

V 1 J 



V^4 = - 



-1 
-1 

V 1 J 



(B5) 



Hence, we can find the solutions by evaluating 

/ Pun \ 



P1010 
P0101 
\ Poooo / 



exp(M 2 2;) 



/ Pun \ 

P1010 
P0101 
\ Poooo / 



(B6) 



z=0 



From the eigenvectors in Eq. ( 1B5I) . we see that the four initial Bell states have zero projection 
to the subspace spanned by (2 and (3 and thus their evolution are not dependent on these 
two parameters. In addition, the elements Piooi^a, ^b, ^>'ai ^'b'j z ) = Pouo( u a, ^b, ^'ai ^'b'j z ) 
and Piioo(^a, ^b, u' a , oj' b \ z) = p* omi {u}A, ^b, u' a , u' b ; z) decay by themselves with the fastest 
decay (4. 

d 



dz 



Pl00l(^A, U>B, w'ai u 'b'i z ) = ~ C4Pl00l(^A, ^' A , Uj' B ] z) 



(B7) 



Putting in the input conditions for the singlet states, i.e. 



1 



Piow(ua,u b ,uj a ,uj b ;0) = po W1 (u a ,ujb,uj' a ,u' b -,0) = -4>(ua,ujb)(P*(uj' a ,u' b ) 
Piooi{va,vb,v a ,u b ;0) = p uo(^A, u B , u' A ,iJ B \ 0) = --<f)(ujA, u) B )<P*(uJ A ,io B ), (B8) 
we obtain the solutions in Eq. fl38|) . 
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APPENDIX C: SOLUTION OF MASTER EQUATION FOR TWO-PHOTON 
STATE IN A COMMON FIBER 



Solving the master equation Eq. 
follows. 



we find that six of the elements are coupled as 
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(C2) 
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with the elements, having a unit of the reciprocal of length, as follows. 



to. 



to 9 



TO4 



TO 

to; 

to' 

m L 

TOn 



-3[/(c A ) 2 + f(u B ) 2 + /K) 2 + /K) 2 ] - 2f{u A )f(u B ) - 2/K)/K) 

+2/(w A )/K) + 2/M/K) + 2f{u B )f{u' A ) + 2f(u A )f(u' B ) 

-3[f(co A ) 2 + /(c B ) 2 + /K) 2 + /(c^) 2 ] + 2f{u A )f{u B ) + 2/K)/(u4) 

+2/(cu A )/K) + 2f(u B )f(u/ B ) - 2f{u B )f{u' A ) - 2f{u A )f{u' B ) 

-3[f{u A f + /K) 2 + /K) 2 + f(io' B ) 2 } + 2f( UA )f( UB ) + 2/K)/K) 

-2/(w A )/K) - 2/(o; B )/K) + 2/(w B )/K) + 2/(o; A )/K) 

4/(w B )/(o/ B ) 

4/K)/K) 

4/(w B )/K) 

Af{u A )f{J B ) 

4/K)/(a; s ) 

4/K)/(u/ B ). (C3) 



To solve the system more conveniently, we introduce the unitary operator, 
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\ _1 _J_ J_ J_ 

\ 2 2^ 

which rotates to the frame where M'(o;o, uo, uo, wo) is diagonal. Applying the transformation 
to M'(uj a ,uj b ,uj a ,uj' b ), 
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where we can see that the 6x6 matrix is expressed as a direct sum of three smaller subma- 
trices, and the matrix elements at u A , wb, uj' a , uj' b are 



(1) 



v 12 

v 22 



(2) 
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-3{[(/( WA ) - f(uj B )} 2 + [(/K) - f(uj' B )} 2 } 
2v / 3[(/(^)-/K)][(/K)-/K)] 

-3[f(u A ) 2 + f(u B ) 2 + f(J A ) 2 + f(u' B ) 2 ] - 2f(u A )f(u B ) - 2f(u' A )f(u' B ) 

+4/( WA )/K) + 4/K)/K) + 4/K)/K) + 4/(w A )/K) 
-3[/K) 2 + f(u B ) 2 + f{u' A ) 2 + f(u' B ) 2 } - 2f(u A )f(cu B ) - 2f(u' A )f(u' B ) 



-2f(u A )f(u' A ) - 2f(uj B )f(uj' B ) - 2f(u B )f(u' A ) - 2f(u A )f(u' 1 
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-3[/( WA ) 2 + f(oo B ) 2 + f(uj' A ) 2 + f(oo' B ) 2 } + 2f(u A )f(u B ) + 2f(u' A )f(u' B ) 
-2f{u A )f(u/ A ) - 2f{u B )f{u' B ) + 2f(u B )f(u' A ) + 2f(u A )f(u' B ) 
A[f(u B )f(u' A ) - f(u A )f(u' B )] 
A[f(u A )f(u B ) - f(u' A )f(u' B )] 

-3[f(co A ) 2 + f(u B f + f{J A ) 2 + /K) 2 ] - 2f{u A )f{u B ) - 2f{J A )f{J B ) 
+2f{u A )f{u' A ) + 2f(u B )f(u' B ) + 2f(u B )f(u' A ) + 2f{u A )f{u' B ) 

4[/K)/K) - fMfW] 

-3[f(u A ) 2 + f{u B f + f{u' A f + /(^) 2 ] + 2f(u A )f(u B ) + 2f(u' A )f(u' B ) 
+2f{u A )f(u' A ) + 2f(u B )f(u/ B ) - 2f{u B )f{u' A ) - 2f{u A )f{u' B ) 



Hence we can solve Eq. (IClh by considering these three submatrices, i.e. 
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If the input is a singlet state as shown in Eq. (150]) . only the first 2x2 submatrix is involved, 
having eigenvalues — & < as follows: 



si — IT 



?2 — — 
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(C8) 



where £2 £1 since V^q \ < 0. We remark that — £1 and —£2 ar e the two eigenvalues 
with the smallest magnitudes in the matrix V. Eq. (1C7I) thus becomes 
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with Vi as presented in Eq. ( 1521) . Then, the output state can be expressed as in Eq. (J5T 
Finally, we point out that or £ 2 to be equal 0, the condition ( V] 

to be satisfied. Considering the function Q(u> A , u>b, uj a , uj' b ) - 



(i) 



(i) 
12 

2 



(1)t/(1) 
11 v 22 



viiVm with a11 



parameters real. In the case when oj' a = u' B , we have Q = only when uja = &b- Otherwise 
when uj' a 7^ uo' B and by fixing uj' a and uo' B we consider dQ/du A = and dQ/du B = 0, giving 
maximum ft at only 2 pairs of real condition, namely u A = u'a an d = w' B , and uja = u' B 
and uj' a = us-, while these two conditions both give Q = 0. 
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